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ABSTRACT 

The unexpected diversity of exoplanets includes a growing number of super-Earth planets, i.e., 
exoplanets with masses of up to several Earth masses and a similar chemical and mineralogical 
composition as Earth. We present a thermal evolution model for a 10 Earth mass planet orbiting 
a star like the Sun. Our model is based on the integrated system approach, which describes 
the photosynthetic biomass production taking into account a variety of climatological, biogeo- 
chemical, and geodynamical processes. This allows us to identify a so-called photosynthesis- 
sustaining habitable zone (pHZ) determined by the limits of biological productivity on the plan- 
etary surface. Our model considers the solar evolution during the main-sequence stage and 
along the Red Giant Branch as described by the most recent solar model. We obtain a large 
set of solutions consistent with the principal possibility of life. The highest likelihood of habit- 
ability is found for "water worlds". Only mass-rich water worlds are able to realize pHZ-type 
habitability beyond the stellar main-sequence on the Red Giant Branch. 
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Introduction 



A pivotal part in the ongoing search for extra- solar planets is the quest to identify p 



a netary 



habitability, i.e., the principal possibility of life. In a previous paper, iKasting et all (| 19931 ) 
presented a one-dimensional climate model to define a zone of habitability (HZ) around the 
Sun and other main-sequence stars that assumed as basic premise an Earth-like model planet 
with a CO2/H2O/N2 atmosphere and that habitability requires the presence of liquid water on 
the planetary surface. 

In the meantime, other definitions of habitable zones have been proposed such as the Galac- 
tic HZ , the UV-HZ, and the photosynthesis-sustaining HZ (pHZ). The Galactic HZ (iLineweaver et al. 



2004)) caters to the requirement that a sufficient amount of heavy elements (notably those con- 



tained in carbon and silicate compounds) must be present fo r the build-up of planets and life, a 
condition easily met in the solar neighborhood. The UV-HZ (IBuccino et a/.i.i2006i : iCuntz et al. 



20081) is based on the premise that no lethal amounts of stellar UV flux is produced (regarding 
life forms assuming carbon-based biochemi stry), a condition l:hat ten ds to favor the environ- 



ment of old main- sequence stars and giants (IGuinan and Ribasl 



appreciable atmospheres, notable significant ozone layers (ISegura et alX 



2002) as we 



2003). 



1 as planets with 



If 

Another definition of habitability first introduced by franck et al.\ pOOOalita) is associated 
with the photosynthetic activity of the planet, which critically depends on the planetary atmo- 
spheric CO2 concentration. This type of habitability is thus strongly influenced by the plane- 
tary geodynamics, encompassing climatological, biogeochemical, and geodynamical processes 
("Integrated System Approach"). This concept has previously been used in studies of fictitious 



plane ts around 47 UMa (,Cuntz et all 120031 : IFranck et alx 120031) and 55 Cnc (|Von Bloh et al 



2003h. as well as detai led studies of observed super-Earth planets in the Gliese 581 system 



von Bloh et a/.Ll2007b|) . The latter investigation showed that Gliese 581c is clearly outside the 



habitable zone, since it is too close to the star, whereas Gliese 58 Id located near the outer edge 
of the ha bitable zone is pro bably habit able, at least for certain types of primitive life forms 



(see also .Selsis et ali l2007h . Moreover, von Bloh et al\ (|2007ah have used this type of model 
to compile a detailed ranking of known star-planet systems regarding the principal possibility 
of life, which by the way led to the conclusion that the Solar System is not the top-tier system 
("Principle of Mediocrity"). 

In case of Ear th-mass planets (1 M m), a detailed investigation of geodynamic habitability 



was presented by iFranck et al.\ (l2000bl) with respec t to the Sun as well as stars of somewhat 



lower and higher mass as central stars. Franck et al\ found that Earth is rendered uninhabitable 
after 6.5 Gyr as a result of plate tectonics, notably the growth of the continental area (enhanced 
loss of atmospheric CO2 by the increased weathering surface) and the dwindling spreading rate 
(diminishing CO2 output from the solid Earth). 

This implies that there is no merit in investigating the future habitability of Earth during the 
post-main- sequence evolution of the Sun, as in the framework of pHZ models, the lifetime of 
habitability is limited by terrestrial geodynamic processes. However, this situation is expected 
to be significantly differ ent for super-E a rth pla nets due to inherent differences compared to 



Earth-mass planets (e.g.. 



Valencia et al. 



2007 a|) . A further motivation for this type of work 



stem s from the ongoing discovery of su per-Earths in the solar neighborhood with the Gliese 



876 ( Rivera et al . 



2005) and Gliese 581 ( Udrvg? al. 



2007|) systems as prime examples. 



In the following, we discuss the definition of the photosynthesis-sustaining habitable zone, 
including the relevant geodynamic assumptions. Next, we describe the most recent model of 
solar evolution that is used as basis for our study. Thereafter, we present our results including 
comparisons to previous work. Finally, we present our summary and conclusions. 



A New Model of Solar Evolution 



A key element of the pre sent study is to consider a s tar akin to the Sun as the central object 



of the star-planet system. 



Schroder and SmithI (|2008|) recently obtained a new model of solar 



evolution that will be adopted in the following. This model is based on a well-tested stellar 
evolution code that allows us to follow the change of solar properties at the main-sequence 
(MS), along the Red Gi ant Branch (RGB) and beyond. It is the Eggleton evolution code in 



the version described by lPols et a/.l(| 199511 1998|) . which has updated opacities and an improved 
equation of state. 

Among other desirable characteristics, the code uses a self-adapting mesh and a proper 
treatment of "overshooting" that has been tested and calibrated with giant and supergiant stars 
in eclipsi ng binary systems. The co de also considers a detailed description of the mass loss 



following ISchroder and Ciinta (|2005|) that has been tested based on a set of well-observed stars 



Schroder and CuntzLl2007|) . Thus it permits an accurate description of the time-dependent solar 



luminosity along the RGB (see Fig. 1). A further consequence of the steadily increasing mass 
loss is the increase of the orbital distances R of any putative planets, given as 7? oc Mq^ with Mq 
as solar mass owing to the conservation of the orbital angular momentum of the planet. 



The solar evolution model by Schroder and SmithI (|2008r ) suggests an age of the Sun as 4.58 
(±0.05) Gyr, and the RGB-tip is reached after 12.167 Gyr, which is also the point in time, 
where our computations are suspended. This model also confirms some well-established facts: 
(1) The MS-Sun has already undergone significant changes, i.e., the present solar luminosity 
exceeds the zero-age value by 0.30 L©, and the zero-age solar radius has been 11% smaller than 
the present value. (2) There was an increase of effective temperature from 5596 K to 5774 (±5) 
K. (3) The present Sun is increasing its average luminosity at a rate of 1% in every 110 million 



years, or 10% over the next bi 



models like the one by Go ug hl (|l98lb and subsequent work 



l ion ye ars. All these findings are consistent with established solar 



During the solar MS, the consequences of evolution for Earth-type planets (as well as other 
types of planets) are extremely slow, compared to natural or human-driven climate changes on 
Earth. Nonetheless, solar-type evolution will force global warming upon any planet, which 
has been the subject of detailed previous i nvestigations both concerning the climatic HZ (e.g.. 



Underwood et al. 



the Su n (e.g.. 



2003 



Jones et al. . 



Franck et al 



2005|) and the photosynthesis-sustain ing habitable zone of 



2000b|) . According to the evolution model by 



Schroder and Smith 



(|2008b . the tip-RGB evolution will be reached with a luminosity of 2730 Lq, an effective tem- 



perature of 2602 K, and a radius of 256 Rq. At that time, the Sun will have lost 0.332 Mq 
of its initial mass. There is an ongoing debate at what point in time a planet originally lo- 
cated at 1 AU, equating 2 1 5 Re-), will be engulfed as a consequence. Contrary to the previous 



mod el by 



Sackmann et al. 



loss. , Schroder and Smith 



(119931) that is based on a less accurate description of the solar mass 



(120081) concluded that such an engulfment will happen during the late 



phase of the solar RGB evolution. In fact, the minimal orbital radius for a planet able to survive 



is found to be about 1.15 AU. 



The evolution of the central star, as well as its effects on planetary orbits, has significant 
consequences for planetary habitability. Thi s property has previo i isly been investig a ted fo r 



different types of climatic habitable zones by iKasting et al\ (|l993|) . lUnderwood et al.\ (|2003h . 



Jones et al\ (l2005h . and others. Furthermore, a previous assessment of the spatial and temporal 



ev olution of clirnatic H Zs for different types of stars beyond the main- sequence has been given 



by 



Lopez et al. 



(|2005|) . They showed that for a 1 Mq star at the first stages of its post-main- 



sequence evolution, the tempor al transit of the HZ is estimated to be several times 10^ yr at 



2 AU and about 10 yr at 9 AU. iLopez et al.\ (|2005b concluded that under these circumstances 



life could develop at distances in the range of 2 to 9 AU in the environment of subgiant or giant 
stars. This view is consistent with our current understanding that terrestrial life existed at least 
as early as 7 x 10^ yr after the Earth formed, which tends to imply that life may be able to form 
over time intervals from 5 x 10^ to 10^ yr. The short-time window 10^ yr) for the origin of 
life is bounded by the last ocean- vaporizing impact and the earliest evidence for life on Earth 
(~3.8 — 3.9x 10^ yr ago). This window rn ight be extended if the origin of life occurred close 



to 3.5 X 10^ yr ago (^Chyba and Hand . 



2005 ) 



The main goal of this study is to investigate habitability in the framework of the photosynthesis- 
sustaining HZ for stars like the Sun with special consideration of the post-main-sequence evo- 
lution. Our study will be focused on super-Earth planets, and we will consider a significant 
set of geodynamic processes. Our findings will also be compared with the previous work by 



Lopez et al. 



( 20051) . 
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Habitability of Super-Earths 



Definition of the photosynthesis-sustaining habitable zone 



To assess the habitability of terrestrial planets, including super-Earth planets, an Earth-system 
model is applied to calculate the evolution of the temperature and atmospheric CO2 concentra- 
tion. On Earth, the carbonate-silicate cycle is the crucial element for a long-term homeostasis 
under increasing solar luminosity. On geological time-scales, the deeper parts of the Earth are 
considerable sink s and sources of carb on. The role of weathering for the Earth's climate was 



first described by IWalker et al\ (|l98l|) . They found that an increase in luminosity leads to a 



higher mean global temperature causing an increase in weathering. As a consequence, more 
CO2 is extracted from the atmosphere, thus weakening the greenhouse effect. Overall the tem- 
perature is lowered and homeostasis is achieved. 

On geological time scales, however, the deeper parts of the Earth are considerable sinks and 
sources for carbo n. As a result, the t ectoni c activity and the continental area change consid- 



erably. Therefore, iTajika and Matsuil (|1992|) have favored the so-called "global carbon cycle". 



In addition to the usual carbonate-silicate geochemical cycle, it also contains the subduction 
of large amounts of carbon into the mantle with descending slabs and the degassing of carbon 
from the mantle at mid-ocean ridges. In particular, the potential of weathering to stabilize the 
surface temperature of a terrestrial planet by a negative feedback mechanism is also strongly 
modulated by the biosphere. 

Our numerical model couples the solar luminosity L, the silicate-rock weathering rate Fwr, 
and the global energy balance to obtain estimates of the partial pressure of atmospheric carbon 
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dioxide /toi' the mean global surface temperature T^uif, and the biological productivity 11 as a 
function of time t (Fig. 2). The main point is the persistent balance between the CO2 (weather- 
ing) sink in the atmosphere-ocean system and the metamorphic (plate-tectonic) sources. This is 
expressed through the dimensionless quantities 

U{t)-fA{t) = /sr(0, (1) 

where /wr(0 = ^wr(0/'^wr.o is the weathering rate, /a(0 =Ac{t)/Acfi is the continental area, 
and /sr(0 = S{t)/So is the areal spreading rate, which are all normalized by their present values 



of Ea rth. Eq. ([T]) can be rearranged by introducing the geophysical forcing ratio GFR (|Volkl . 



19871) as 



U{TsM,Pco2) = ^ =: GFR{t). (2) 

JA 

Here we assume that the weathering rate depends only on the global surface temperature and the 
atmospheric CO2 concentration. For the investigation of a super-Earth under external forcing, 
we adopt a model planet with a prescribed continental area. The fraction of continental area 
relative to the total planetary surface c is varied between 0. 1 and 0.9. 

The connection between the stellar parameters and the planetary climate can be formulated 
by using a radiation balance equation 

^[l-a(r,urf,/'co2)] = 4//j(rsurf,/'co2), (3) 
where L denotes the stellar luminosity, R the pla netary distance, a the planetary albedo, and Ir 



the outgoing infrared flux of the planet. Following iWilliamsi(|1998|) Ir has been approximated by 



a third order polynomial and a by a second order polynomial. These approximations have been 
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deri ved from 24, OO P runs of a radiation-convection model by iKasting and AckermanI (|l986h 



and 



KastingI (Il988b . They are valid in a range of 10 ^ bar < /'c02 < 10 bar. The Eqs. © 



and dl]) constitute a set of two coupled equations with two unknowns, T^uif and Pco2, if the 
parameterization of the weathering rate, the luminosity, the distance to the central star and the 
geophysical forcing ratio are specified. Therefore, a numerical solution can be attained in a 
straightforward manner. 

The photosynthesis-sustaining HZ is defined as the spatial domain of all distances R from 
the central star, e.g., the Sun, where the biological productivity is greater than zero, i.e.. 



pHZ := {R\n{Pco2{R.t),TsuAR,t))>0}. 



(4) 



In our model, biological productivity is considered to be solely a function of the surface tem- 
perature and the CO2 partial pressure in the atmosphere. Our parameterization yields maxi- 
mum productivity at T^urf = 50 ° C and zero productivity for Tsurf < 0°C or Tjurf > 100°C or 



Pc02 < 10 ^ bar (|Franck et ali 



2000al) . A photosynthesis-based biosphere of a super-Earth 



may, however, use methane to produce CO2, because hydrogen is less likely to escape to space. 
The inner and outer boundaries of the pHZ do not depend on the detailed parameterization of 
the biological productivity within the temperature and pressure tolerance window. Hyperther- 
mophilic life forms can tolerate temperatures somewhat above 100°C. However, these chemoau- 



totrophic organisms are outside the scope of this study. 



12 



Silicate Rock Weathering 

Weathering plays an important role in Earth's climate because it provides the main sink for 
atmospheric carbon dioxide. The overall chemical reactions for the weathering process are 

C02+CaSi03 ^ CaC03 + Si02, 
C02+MgSi03 ^ MgC03 + Si02. 

The total process of weathering embraces (1) the reaction of silicate minerals with carbon diox- 
ide, (2) the transport of weathering products, and (3) the deposition of carbonate minerals in the 
oceanic crust. The available thickness of crust where CaC03 is stable in the presence of silicate 
scales inversely with the thermal gradient and hence inversely with surface gravity. Therefore, 
there may be a problem for storing carbonates in the crust of super-Earth planets. Additionally, 
there is an exchange with the mantle via alteration of the oceanic crust. 

When combining all these effects, the normalized global mean weathering rate /wr can be 
calculated as 

^H+ \ ^'^ f T'surf — 7surf,0 



following 



Walker et al. 



(|1981|) . Here the first factor reflects the role of the CO2 concentration 
in the soil, Psoib with aH+ as the activity of H"*" in fresh soil- water that depends on P^oii and 
the global mean surface temperature T^urf- The quantities 0^+ ^i^^ ^surf.o are the present-day 
values for the H+ activity and the surface temperature, respectively. The activity aii+ is itself a 
function of the temperature and the CO2 concentration of the soil. The concentration of CO2 in 
the soil water [C02(aq)] can be obtained from the partial pressure of CO2 in the soil according 
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to 



[C02{aq)] = KuPsoil 



(6) 



where is Henry's law constant. We assume that [C02(aq)] = [H2CO3]. H2CO3 dissociates 
in two steps which are 



H2CO* 



HCO, 



H+ + HCO; 



H+ + CO 



2- 



The corresponding concentrations can be calculated from the law of masses as 



fHCO. 



CO' 



[H+. 
K1K2 

WW 



KuPsoih 



soil; 



(7) 
(8) 



where ^1 and K2 are (temperature dependent) equilibrium constants. An additional constraint 
for the concentrations is given by the charge balance 



+1 



[HC03]+2[C02-] + [OH 



(9) 



Here [H+] and therefore can be derived from a simultaneous solution of Eqs. (O to dH) as 
a function of Pgoii- The sulfur content of the soil can be taken into account analogously. The 
equilibriu m constants for the chemical activities of the carbon and sulfur systems involved are 



taken from lStumm and MorganI (Il98lh . Note that the sulfur content of the soil also contributes 
to the global weathering rate, but its influence does not depend on the temperature. It can be 
regarded as an overall weathering attribute that has to be taken into account for the estimation 
of the present-day value. 



14 



For any given weathering rate, the surface temperature and the CO2 concentration of the 
soil can be calculated in a self-consistent manner. A nii is assumed to be linearly related to the 



terrestrial biological productivity U (see IVolkl Il987h and the atmospheric CO2 concentration 



Pco2- Thus we find 



■fsoil 



^CO2,0^ _^ Pco2 



(10) 



^soil,0 \ Aoil.O / Aoil,0 

where Psoii.o^ no and PCO2.0 are again present-day values. Note that the present-day concentra- 
tion of CO2 in the soil is ten times the present-day concentration of CO2 in the atmosphere, i.e., 

^soil =10 /'soiLO- 



Thermal Evolution Model 



Parameterized convection models are the simplest models for investigating the thermal evo- 



lution of terrestrial planets and satellites. They have s uccessfully been 



lution of Mercury, Venus, E arth, Mars, and the Moon (jStevenson et al. 



applied to the evo 



1983 



Sleepl . 



2000 ) 



Franck and Bounamal (| 19951) have investigated the thermal and volatile history of Earth and 



Venus in the framework of comparative planetology. The internal st ructure of mass i ve teir es 



trial planets with one to ten Earth masses has been investigated by [Valencia et al\ (l2006h to 



obtain scaling laws for the total radius, mantle thickness, core size, and average density as a 
function of mass. Further scaling laws were found for different compositions. We will use 
such scaling laws for mass-dependent properti es of our 10 super-Earth m odel as well as for 
mass-independent material properties given by iFranck and Bounamal (Il995h (see Tab.[T]). 



The thermal history and future of a super-Earth has to be determined to calculate the spread- 
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ing rate for solving Eq. ©. A parameterized model of whole mantle convection inc 



the volatile ex change between the mantle and surface reservoirs (|Franck and Bounama . 



uding 



1995 



Franckl . ll998|) is applied. Assuming conservation of energy, the average mantle temperature r„, 



can be obtained by solving 

^-npc{Rl-Rl)^ = -AnRlq^ + '^-'KE{t){Rl-Rl), (11) 

where p is the density, c is the specific heat at constant pressure, qm is the heat flow from the 
mantle, E{t) is the energy production rate by decay of radiogenic heat sources in the mantle 
per unit volume, and Rm and Rc are the outer and inner radii of the mantle, respectively. The 
radiogenic heat source per unit volume is parameterized as 

E{t) = Eoe-^ (12) 

where X is the decay constant and the constant Eq is obtained from the present heat flux of 
qm = 0.07 Wm^^ for an Earth-size planet at 4.6 Gyr. 

The mantle heat flow is parameterized in terms of the Rayleigh number Ra as 

k{T,n-T,M) f V 

1m = (13) 



Rm ~Rc V R^crit / 

with 

^ ga{T,n-T,^rf){Rm-Rc? 

KV ' 

where k is the thermal conductivity, Racrit is the critical value of Ra for the onset of convection, 
(3 is an empirical constant, g is the gravitational acceleration, a is the coefficient of thermal 
expansion, K is the thermal diffusivity, and v is the water-dependent kinematic viscosity. The 
viscosity v can be calculated with the help of a water fugacity-dependent mantle creep rate. It 

16 



strongly depends on the evolution of the mass of mantle \ yater Mw, and the man t 



T,n, i.e., V = y{T,n,Mw) and is parameterized according to 



Franck and Bounamal (Il995h . 



^temperature 



The evolution of the mantle water can be described by a balance equation between the 
regassing flux Freg and outgassing flux Fout as 



dt 



— ^reg ^out 



/basPbas^i?bas5''^H20 4 



dmfw^i 



(15) 



where /bas is the water content in the basalt layer, Pbas is the average density, Jbas is the av- 
erage thickness of the basalt layer before subduction, S is the areal spreading rate, dm is the 
melt generation depth and is the outgassing fraction of water. i?H20 is the regassing ratio of 
water, i.e., the fraction of subducting water that actually enters the deep mantle. The average 
thickness of the basalt layer as well as the melt generation depth scale inversely with surface 
gravity g, i.e., Jbas °^^/g and dm°^ 1/^. The pressure closing of cracks in the deeper parts of 
the basalt layer scales also inversely with g and thus reduces the storage capacity of volatiles for 
a super-Earth planet. Therefore, the ratio Freg/Fout is independent of g. According to Eq. (fTS]) 
gravity influences only the time scale of mantle water evolution. Therefore, as a first approxi- 



mation the melt generation depth d m does not depend on mantle 



is a temperature dependence of dm (|McKenzie and Bickld . 



1988 



temperature. However , there 



Langmuir et al. 



1992 ). The 



regassing ratio depends linearly on the mean mantle temperature Tm that is derived from the 
thermal evolution model via 



Rn^oiTm) = RT-{Tm{0)-Tm)+Rn2O,Q. 



(16) 



The factor Rj is adjusted to obtain the correct modem amount of surface water (one ocean 
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mass) for an Earth-size planet and 7?h2O,0 is fixed at 0.001. This value is obviously very low 
at the beginning of the planetary evolution because of the enhanced loss of volatiles resulting 
from back-arc volcanism at higher temperatures. 

The areal spreading rate 5 is a function of the average mantle temperature Tm, the sur- 
face temperature T^ur^,^ the h eat flow from the mantle qm, and the area of ocean basins Aq 
jTurcotte and Schubert!, 1982,) . given as 

S ^ ^^''^^0 (17) 

( 7m T'surf ) ^ 

In order to calculate the spreading rates for a planet with several Earth masses, the planetary 
parameters have to be adjusted accordingly. We assume 

0.27 





and with Rp as planetary radius, see dValencia et al. . 



(18) 



20061). The total radius, mantle thickness. 



core size and average density are all functions of mass, with subscript © denoting Earth values. 
The exp onent of 0.27 has been obtained for super-Earths (M > IM®), and has already been 



used by 



von Bloh et al.\ (|2007b|) in their models of Gliese 581c and 58 Id. The values of R^, Rc, 



Aq, the density of the planet, and the other planetary properties are also scaled accordingly. 

The source of CO2 to the atmosphere is expressed in mass of carbon outgassed at the spread- 
ing zones, Csr °^ S. It has to be converted to an equivalent concentration of CO2 in the atmo- 
sphere. This can be done by the following equation 

^CO, = 7^^C,r, (19) 

where ^ucoa ^i^d /jq are the molar weights of CO2 and C, respectively. The mass dependent pre- 
factor g/R^ scales as M^'^-^^ ^ M*^ and has therefore been neglected in our study. Therefore 
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the conversion does not depend on the planetary mass and the spreading rates S can be directly 
used to calculate /si- in Eq. ([T])- 

In Tab. \T\ we give a summary of the selected values for the parameters used in the thermal 
evolution model of the 10 super-Earth planet, while also depicting an Earth-size planet for 



comparison. According to [Valencia et all (|2007bh . we assume that a more massive planet is 



likely to convect in a plate tectonic regime similar to Earth. Thus, the more massive the planet 
is, the higher the Rayleigh number that controls convection, the thinner the top boundary layer 
(lithosphere), and the faster the convective velocities. This is the so-called boundary-layer limit 
of convection. From this limit it follows that the interior of a super-Earth is always hotter and 
less viscous than that of a Earth-mass planet. Nevertheless, friction is the rate-limiting process 
for subduction. Increa sing the planetary radius acts to decrease the ratio between driving forces 



and resistive strength (|0'Neill and LenardicL l2007h . Thus a super-sized Earth might be in an 



episodic or stagnant lid regime. 

In a first order approximation, we assume a fixed thickness of the basalt layer and melting 
depth corresponding to relatively low values. Furthermore, the initial amount of water ^^,(0) 
scales linearly with the planetary mass. However, his might be an underestimate because more 
massive planets tend to accrete more volatiles. 
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Results 



Habitability based on the integrated system approach 



approach that has p 


1 

reviously been used in various other planetary studies (e.e..iFrancke?a/.. 


2000b: 


Cuntz et al. 


2003; 


von Bloh et al. 


2003; 


Franck et al. . 


2003; 


von Bloh et al. 


2007a 


b). 



The simulations have been carried out for a 10 M® mass super-Earth with a fixed relative conti- 
nental area c varied from 0.1 to 0.9. Fig. 3 shows the behavior of the photosynthesis-sustaining 
habitable zone (pHZ) of the Sun for a 10 M® super-Earth planet. The age domain beyond 
1 1 Gyr that also includes the post-main- sequence evolution is depicted in Fig. 4. The width 
of the pHZ during the main-sequence evolution is found to be approximately constant, but for 
higher ages, it increases over time and moves outward, a phenomenon most noticeable beyond 
11.5 Gyr. For example, for ages of 11.0, 11.5, 12.0, and 12.1 Gyr, the pHZ is found to extend 
from 1.41 to 2.60, 1.58 to 2.60, 4.03 to 6.03, and 6.35 to 9.35 AU, respectively. 

At relatively high ages, habitable solutions are identified as water worlds, if the Sun as 
central star has reached the RGB. The reason is that planets with a considerable continental area 
have higher weathering rates that provide the main sink of atmospheric CO2. Therefore, such 
planets are unable to build up C02-rich atmospheres which prevent the planet from freezing or 



allowing photosynthesis-based life. This result is co nsistent with previous findings for ^ Earth 



mass planets around the Sun or stars of similar masses (iFranck et al. 



2000bl : 



Cuntz al.u2003 



Note that the partial pressure of carbon dioxide in the planetary atmosphere is determined 
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by the equilibrium of sources and sinks. The sources are given by volcanic outgassing, while 



the sinks are given by the weat hering of si 



studies of 1 planets (e.g., 



Franck et al. 



i cates o n the continents. As previously found in 



2000b|) . the rate of outgassing is monotonously 



decreasing with age because of the decay of long-lived radionuclides and the loss of the initially 
available accretion energy at the planetary surface. This process starts just after completion 
of the planetary accretion both for an initially habitable and uninhabitable planet. A planet 
beyond the outer edge of the pHZ is completely frozen and thus no weathering will occur on 
the continents. Therefore, all CO2 is accumulated in the atmosphere. If the planet becomes 
habitable due to the increase of the luminosity of the central star, weathering starts and a new 
equilibrium of atmospheric CO2 is established as a consequence. 

Furthermore, the interior of a planet with a relatively low mass is known to cool down 
more rapidly. Therefore, such a planet initially beyond the outer edge of the habitable zone 
will not become habitable at a later stage because of the failure to provide a sufficiently dense 
atmosphere. In contrast, a super-Earth planet might become habitable, depending on the relative 

i f snow ball planets as a possible 



size of the continental area. In a recent study, the impor tance o 
source of water-rich terrestrial planets was elucidated by iTajikal (|2008h . although the main focus 
of this paper was the assessment of internal oceans. 

Super-Earth-type water worlds are even able to realize pHZ-type habitability beyond solar- 
type main-sequence evolution. Any model where mantle vents its water will end up as a water 
world super-Earth. The height of ridges, volcanos, and mountains scale with lithosphere thick- 
ness and hence with l/g. As the central star evolves, its pHZ expands outward, and moves 
further away from the star, particularly for stellar ages beyond 1 1.8 Gyr (see Fig. 4). Similar to 
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the climatic HZ ("see .Lopez et aL for details), the pHZ acts like a shell that sweeps progressively 
outward over a wide range of distances from the star. This results in a significant decrease of the 
duration of the transit of the habitable zone for any planet located beyond 1.5 AU (see Fig. 5). 
We find that for water worlds with c = 0.1, the duration of the transit of the pHZ at 2, 3, and 
5 AU is given as 3.7, 0.25, and 0.10 Gyr, respectively, whereas for planets at 10 and 20 AU, 
much smaller durations of the transit are identified. 

Figs. 3 and 4 also depict various orbital distances of planets originally located between 1 AU 
and 5 AU. Note that the se orbital distances do not change during the stellar main-sequence 



stage, i.e., below 10 Gyr (|Schroder and SmithLl2008h . owing to the lack of significant mass loss 



and in the absence of significant planet-planet interaction as typically encountered in multiple 
planetary systems. Thereafter, the orbital distances R of any planet increases following R Mq^ 
with Mq as mass owing to the conservation of planetary orbital angular momentum. 

Comparison with previous results 

The existence of habitability around stars that h ave evolved away f rom the main-sequence has 



already been the topic of previous inv estigations. 



Lopez et al. 



climatic HZ 



jased on the conce pt of iKasting et al.\ (| 1993b for stars of different masses, includ 



(|2005r) studied the behavior of the 



ing the Sun. iLopez et all (120051) assume a HZ based on conservative limits of habitability. The 



inner limit of their conservative estimate is set by the lowest temperature at which the liquid- 
solid phase change of water occurs. The estimat e of the outer l imit a ssumes the existence of 



a greenhouse effect involving CO2 and H2O gas (IKasting et al. 



1993 ). The less conservative 



definition extends the outer edge of the limit of habitability to as large as 2.4 AU, largely de- 



22 



and 



Mischna et al. 



pen ding on the rad i ative p roperties of the CO2 ice clouds; see forget and PierrehumbertI (|1997f ) 
( 2000 ) for detailed studies. 



Akin to the pHZ previously discussed. iLopez et al\ (|2005[) found that for the Sun during its 



evolution the climatic HZ acts like a shell that sweeps progressively outward over a wide range 
of distances from the star. The duration of the transit during which the HZ passes over a planet 
located at 1 AU from the star was found to be on the order of 10^ yr. After the star leaves the 
main- sequence, the climatic HZ progressively moves to 2 AU. The duration of the transit at 
this location is approximately 10^ yr. A plateau is observed in the curve up to 9 AU (for the 
conservative limits) and up to 13 AU (for the less conservative limits), where the durations of 
habitable conditions lasts from a few to several times 10^ yr. At 10 AU, the duration is smaller, 
about 10^ yr. At 15 AU from the star the duration of habitable conditions lasts more than 10^ yr. 



and at the largest distances considered in the study by lLopez et a/.l(|2005h the duration gradually 



decreases. 



Note that the model of solar evolution considered in the lLopez et al\ study is that by lMaeder and Meynet 



by 



(119881). Nonetheles s , their results would be quite simil ar if they had used the subse quent model 



Sackmann et ah (il993|) or the very recent model by iSchroder and Smith 



(20081). The reason 



is that the outcome of the lLopez et al\ study is much more dependent on the choices made con- 
cerning the upper and lower limits of the climatic HZ, mostly connected to the treatment of the 
CO2 atmospheres (i.e., radiative properties, cloud coverage, etc.) than to the adopted model of 
solar evolution. 



I 

Fig. 6 shows the comparison between the work by iLopez et al. 



(2005) and our current re- 



sults. We find that for water worlds (c = 0.1), the transit times for photosynthesis-sustaining 
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habitability (pHZ) for planets at a given reference distance from the star is re 



lar to the re sults obtained for the conservative climatic HZ (IKasting et al. 



19931) 



atively simi- 
adopted by 



Lopez et al.l albeit the transit times in our study are typically lower by a factor of up to 1.5. For 



example, the durations of the transit concerning pHZ-type habitability for water worlds at 2, 3, 
and 5 AU are identified as 3.7, 0.25, and 0.10 Gyr, respectively, whereas for planets at 10 and 
20 AU, the durations of the transit found are as low as 27 and 9 Myr, respectively. 

However, especially at distances beyond 2 AU, significantly smaller transit times are en- 
countered for planets with larger continental areas in terms of all stellar distances, a result con- 
sistent with previous findings. For c = 0.5, the transit time of the pHZ drops beneath 1 Gyr for 
planets located at 1.8 AU. For planets with a relative continental area of c = 0.9, also referred 
to as "land worlds", no significant photosynthesis-sustaining habitability is found for planets 



beyond 1.5 AU. 



Summary and Conclusions 



' 1 

has previously been used in various theoretical planetary studies ( e.s... Franck et aU 


1 

,2000b; 


Cuntz et al. . 


2003; 


von Bloh et al. . 


2003; 


Franck et al. . 


2003; 


von Bloh et al. . 


2007a t 


)). This 



work is motivated by the quest to identify habitability outside the Solar System as well as the 



ongoing discovery of s uper-Earths in the solar neighborhood with the Gliese 876 (|Rivera et al. 



20051) and Gliese 581 (lUdrv al. . 



2007|) systems as prime examples. 



In agreement with previous studies, it is found that photosynthesis-sustaining habitability 
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strongly depends on the planetary characteristics. For planets of a distinct size, the most impor- 
tant factor is the relative continental area. Habitability was found most likely for water worlds, 
i.e., planets with a relatively small continental area. For planets at a distinct distance from the 
central star, we identified maximum durations of the transit of the pHZ. A comparison of planets 
with different masses revealed that the maximum duration of the transit increases with planetary 
mass. Therefore, the upper limit for the duration of the transit for any kind of Earth-type planet 
is found for most massive super- Earth planets, i.e., 10 M®, rather than 1 Mm planets, wh ich are 



rendered uninhabitable after 6.5 Gyr, as previously pointed out by 



Franck et al.\ (l2000b ) 



Our study forwards a thermal evolution model for a 10 M® super-Earth orbiting a star 



akin to the Sun. The ca 



culations consider updated models of solar evo lution obtained by 



Schroder and SmithI (120081) with a detailed mass loss description provided by iSchroder and Cuntz 



(|2005r ). The latter is relevant for the change of luminosity along the Red Giant Branch as well 
as the increase of the orbital distances of any putative planets during that phase. By employing 
the integrated system approach, we were able to identify the sources and sinks of atmospheric 
carbon dioxide on the planet, allowing us to describe the photosynthesis-sustaining habitable 
zone (pHZ) determined by the limits of biological productivity on the planetary surface. 
Concerning the pHZ, we identified the following properties: 

(1) Geodynamic solutions are identified for different solar ages, including the RGB phase. The 
pHZ increases in width over time and moves outward. For example, for ages of 1 1 .0, 11 .5, 12.0, 
and 12.1 Gyr, the pHZ is found to extend from 1.41 to 2.60, 1.58 to 2.60, 4.03 to 6.03, and 6.35 
to 9.35 AU, respectively. 

(2) Habitable solutions at large ages, especially for the subgiant and giant phase, are water 
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worlds. This also means that the possibility of water worlds in principle results in an extension 
of the outer edge of habitability. The reason is that planets with a considerable continental area 
have higher weathering rates that provide the main sink of atmospheric CO2. Therefore, such 
planets, contrary to water worlds, are unable to build up COi-rich atmospheres that prevent the 
planet from freezing or allowing photosynthesis-based life. 



(3) Th e total duration of the transit of the habitable zone is similar to the pr edictions bvlLopez et al. 



(|2005r) based on the conservative limits of the climatic HZ obtained by 



Kasting et al. 



(119931) 



For water worlds with c = 0.1, the transit times of the pHZ at 2, 3, and 5 AU are obtained as 
3.7, 0.25, and 0.10 Gyr, respectively, whereas for planets at 10 and 20 AU, much smaller transit 
times are found. 



Our results are a further motivation to consider super-Earth planets in upcoming or proposed 
planet search missions such as Kepler, TPF or Darwin. Moreover, our results can also be viewed 
as a reminder to seriously contemplate the po ssibility of habitable planets around red giants, as 



previously pointed out by iLopez et al.\ (|2005|) and others. For central stars with a higher mass 



than the Sun, a more rapid evolution will occur that will also affect the temporal and spatial 
constraints on planetary habitability when the central stars have reached the RGB. 
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Table 1 : Parameter values for the evolution model of the mantle temperature and water 



Parameter 


Value 


Unit 


Description 




IMe 


10 Me 








5x 10^ 


5x 10^ 


m 


Average thickness of the basalt layer 


/bas 


0.03 


0.03 




Mass fraction of water in the basalt layer 


Pbas 


2,950 


4,569 


kg m-^ 


Density of the basalt 


fw 


0.194 


0.194 




Degassing fraction of water 


dm 


40 X 10^ 


40 X 10^ 


m 


melting depth 


k 


4.2 


4.2 


J s 1 ml K 1 


Thermal conductivity 


Rc 


3,471 X 10^ 


6,463 X 103 


m 


Inner radius of the mantle 


Rm 


6,271 X 10^ 


11,667 X 10^ 


m 


Outer radius of the mantle 




4.2 X 10^1 


14.04 X 10^2 


kg 


Initial amount of mantle water 




3,000 


3,000 


K 


Initial mantle temperature 




3.36 X 10-4 


3.36 X 10-4 


mol J-i 


Henry's law constant at 25 °C 




-6.3 


-6.3 




Equilibrium constant in Eq. (|7]) at 25°C 


l0gi^2 


-10.3 


-10.3 




Equilibrium constant in Eq. ([8]) at 25°C 


K 


10-*^ 


10^ 


2 -1 

m s 


Thermal diffusivity 


pc 


4.2 X 10^ 


4.2 X 10^ 


Jm 3 K 1 


Density x specific heat 


Rt 


2.98 X 10-4 


2.98 X 10-4 


K 1 


Temperature dependence of regassing ratio 


a 


3 X 10-5 


3 X 10-5 


K 1 


Coefficient of thermal expansion 


P 


0.3 


0.3 




Empirical constant in Eq. (fT3l) 


R^crit 


1,100 


1,100 




Critical Rayleigh number 


X 


0.34 


0.34 


Gyr 1 


Decay constant 


Eo 


1.46 X 10-^ 


1.46 X 10-^ 


J s 1 m^ 


Initial heat generation per time and volume 


8 


9.81 


28.26 


-2 

m s 


Gravitational acceleration 



34 




-1 







4 6 8 

Age (Gyr) 



1 



1 2 



-0 



Figure 1 : Solar evolution model following 



Schroder and Smithl(|2008|) . depicting the luminosity 



(solid line), the effective temperature (dashed line), and the mass (dash-dotted line). 



35 




Figure 2: Box model of the integrated system approach adopted in our model. The arrows 
indicate the different types of forcing, which are the main feedback loop for stabilizing the 
climate (thick solid arrows), the feedback loop within the system (thin solid arrows), and the 
external and internal forcings (dashed arrows). 
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Figure 3: The photosynthesis-sustaining habitable zone (pHZ) of the Sun for a super-Earth 
planet also considering the solar evolution along the Red Giant Branch. The colored areas 
indicate the extent of the pHZ for different relative continental areas. The solid line depicts the 
orbital distance of a planet originally located at 1.0 AU. 
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Figure 4: Same as Fig.[3l but now zoomed-in at the Red Giant Branch evolution beyond 1 1 Gyr. 
We also show the distance evolution of various fictitious planets, originally located at 1 .0 AU 
(solid line), 1.5 AU (dotted line), 2.0 AU (dash-dotted line), 3.0 AU (dashed line), and 5.0 AU 
(dash-double dotted line), respectively. 
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Figure 5: Duration of the transit of the pHZ for a super-Earth as function of the distance from 
the central star R for different relative continental areas. Model (a) only includes the stellar 
main-sequence evolution (10 Gyr), whereas Model (b) also includes the evolution along the 
Red Giant Branch . The g ray area indicates the result based on the geostatic approximation 



(e.g., iFranck et al. 



2000b|) that is mostly enveloping the colored area. The dash-dot t ed lin e 



indicates transit times of the general habitable zone given by the model of 



Kasting et al. 



(119931) 



39 



10' 



10' 



o 10 

Q. 
(/) 

CD 

^ 10 



10 



-1 



-2 



-3 







10 15 

R (AU) 



20 



25 



Figure 6: Duration of the transit of the pHZ for super-Earth planets as function of the plane- 
tary distance R from the center of the star. We depict the results for "water worlds" (c = 0.1; 



dashed line) and " 



Lopez et al. 



Kasting et al. 



and worlds" (c = 0.9; dotted line). In addition, we also show the results by 



( 2005h for terrestrial planets (solid lines) based on the concept of habitability by 



19931) . Here the lower line refers to the duration of habitability based on their 
conservative definition of the HZ, whereas for the upper line the less conservative definition has 
been adopted for the outer limit of the HZ, while the inner limit of the HZ was left unchanged. 
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